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(54) Abstract Title: Investigation of subterranean reservoirs 

(57) A method of investigating subterranean strata (14), particularly in the seabed, comprises the steps of 

deploying an electric dipole transmitter antenna (16); deploying an electric dipole receiver antenna (15) at 
a predetermined offset from the transmitter; applying an electromagnetic wave field to the strata and 
detecting the wave field response; extracting phase information from the response; repeating the 
procedure with the transmitter and/or receiver in different locations for a plurality of transmissions; and 
using the phase information from the plurality of transmission responses to determine the presence 
and/or nature of any reservoir (18) present; in which method the transmission wavelength X through the 
overburden is between 0.1s and 10s where s is the distance from the seabed to the reservoir, the offset 
between transmitter and receiver is between 0.5^ and 10X, and the transmission frequency is from 0.01 Hz 
to 1kHz. A refracted wave response provides useful information in the wave field response. 
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2413187 



Method and Apparatus f or Determining 
the Nature of Subterran ean Reservnirg 

The present invents relates to a method and apparatus for determining the nature of 
submarine and subterranean reservoirs. The invention is particularly suitable for 
determining whether a reservoir, whose approximate geometry and location are 
known, contains hydrocarbons or water, and also for detecting reservoirs with 
particular characteristics. 

Currently, the most widely used techniques for geological surveying, particularly in 
sub-marine situations, are seismic methods. These seismic techniques are capable of 
revealing the structure of the subterranean strata with some accuracy. However 
whereas a seismic survey can reveal the location and shape of a potential reservoir, it 
can normally not reveal the nature of the reservoir. 

The solution therefore is to drill a borehole into the reservoir. However, the costs 
involved in drilling an exploration well tend to be in the region of £25m and since the 
success rate is generally about 1 m 10, this tends to be a very costly exercise. 

It is therefore an object of the invention to provide a system for locating a 
subterranean reservoir and for determining, its nature with greater certainty, without 
the need to sink a borehole. 



It has been appreciated by the present applicants that while the seismic properties of 
hydrocarbon filled strata and water-filled strata do not differ significantly their 
electromagnetic resistivities do differ. Thus, by using an electromagnetic surveying 
method, these differences can be exploited and the success rate in predicting the nature 
of a reservoir can be increased significantly. This represents potentially an enormous 

l~i C- f OAT r , 



cost saving. 



Consequently, a method and apparatus embodying these principles from the basis of 
the present applicant's International British patent application No. WO 01/57555. 

This contemplates a method of determining the natnre of a subterranean reservoir 
whose approximate geometry and location are known, which comprise.,: applying a 
..me varying electromagnetic field ,„ ,he strata containing the reservoir; detecting the 
electromagnetic wave field response; seetang in «he wave field response, a component 
representing a refracted wave from the hydrocarbon layer; and determining the content 
of the reservoir, based on the presence or absence of a wave component refracted by 
the hydrocarbon layer. 

It also contemplates a method for searching for a hydrocarbon containing subterranean 
reservon- which comprises: applying a time varying electromagnetic field to 
subterranean strata; detecting the electromagnetic wave field response; seeking in the 
wave field response, a component representing a refracted wave; and determining the 
presence and/or nature of any reservoir identified based on the presence or absence of 
a wave component refracted by hydrocarbon layer. 

It further contemplates an apparatus for determining the nature of a subterranean 
reservoir whose approximate geometry and location are known, or for searching for a 
hydrocarbon containing subterranean reservoir, the apparatus comprising: means for 
applying a time varying electromagnetic field to the strata containing the reservoir- 
means for detecting the electromagnet, wave field response; and means for seeking' 
in the wave field response, a component representing a refracted wave, thereby 
enabling the presence and/or nature of a reservoir to be determined. 

A refracted wave behaves differently, depending on the nature of the stratum in which 
it is propagated. In particular, the propagation losses in hydrocarbon stratum are much 
lower than in a water-bearing stratum while the speed of propagation is much higher 
Fhus, when an oil-bearing reservoir is present, and an EM field is applied, a strong and 



rapidly propagated refracted wave can be detected. This may therefore indicate the 
presence of the reservorr or its nature if its presence is already known. 

Electromagnetic surveying techniques in themselves are known. However, they are 
not w ]d ely used in practrce. In general, the reservoirs of interest are about 1 km or 
more below the seabed. In order to carry out electromagnetic surveying as a stand 
alone technique in these conditions, with any reasonable degree of resolution, short 
wavelengths are necessary. Unfortunately, such short wavelengths suffer from very 
h.gh attenuation. Long wavelengths do not provide adequate resolution. For these 
reasons, seismic techniques are preferred. 



lues cannot 



However, while longer wavelengths applied by electromagnetic techniq, 
provide sufficient information to provide an accurate indication of the boundaries of 
the various strata, if the geological structure is already lenown, they can be used to 
determine the nature of a particular rdentified formation, if the possibilities for the 
nature of that formation have significantly differing electromagnetic characteristics 
The resolutron is not particularly important and so longer wavelengths which do not 
suffer from excessive attenuation can be employed. 

The resistivity of seawater is about 0.3 ohm-m and that of the overburden beneath the 
seabed would typically be from 0.3 to 4 ohm-m, for example about 2 ohm-m 
However, the resistivrty of a hydrocarbon reservorr is likely to be about 20-300 ohm- 
m. Typrcally, therefore, the resistivity of a hydrocarbon-bearing formation will be 20 
to 300 tunes greater than that of a water-bearing formation. This large difference can 
be explo.ted using the techniques of the present invention. 

The electrical resistivity of a hydrocarbon reservoir normally is far higher than the 
surrounding material (overburden). EM-waves attenuate more rapidly, and travel 
slower inside a low resistivity medrum, compared to a high resistivity medium 
Consequently, hydrocarbon reservoir will attenuate EM-waves less, compared to a 



lower resistivity overburden. Furthermore, the EM-wave speed will be higher inside 
the reservoir. 

Thus, an electric dipole transmitter antenna on or close to the sea floor induces 
electromagnet. EM fields and currents in the sea water and in the subsurface strata 
hi the sea water, the EM-fields are strongly attenuated due to the high conductivity in 
the sahne environment, whereas the subsurface strata with less conductivity causes 
less attenuation. If the frequency is low enough (in the order of 1 Hz), the EM energy 
is able to penetrate deep into the subsurface, and deeply buried geological layers 
havmg higher electrical resistivity than the overburden (as e.g. a hydrocarbon filled 
reservoir) will affect the EM-wave, Depending on the angle of incidence and state of 
polarisation, an EM wave incident upon a high resistive layer may excite a ducted 
(guided) wave mode in the layer. The ducted mode is propagated laterally along the 
layer and leaks energy back to the overburden and receivers positioned on the sea 
floor. In the present application, such a wave mode is referred to as a "refracted 



wave". 



The distance between the EM source and a receiver is referred to as the offset. Due to 
the fact that a refracted wave in a hydrocarbon-bearing formation will be less 
attenuated than a direct wave in seawater (or in the overburden), for any given H/C 
bearing formation, there will be a critical offset at which the refracted wave and the 
direct wave will have the same signal strength. This may typically be about two to 
three times greater than the shortest distance from the source (or receiver to the H/C 
bearing formation). Thus, when the offset is greater than the critical offset, the radial 
EM waves that are refracted into, and guided through the reservoir, will pay a major 
contribution to the received signal. The receiver signal will be of greater magnitude 
and arrive earlier (i.e. have smaller phase) compared to the case where there is no HC 
reservoir. In many cases, the phase change and/or magnitude change recorded at 
distances greater than the critical offset, may be directly used for calculation of the 
reservoir resistivity. Furthermore, the reservoir depth may be inferred from the critical 



offset and/or the phase and magnitude slopes for 



various source - receiver offsets. 



According to one aspect of the present invention, there is piovided a method of 
tnvesfgating snbleiranean slrala which comprises: deploying an eieclric dipo.e 
transmitter antenna; defying an eleclric dipole receiver anlenna a, a predefined 
offset d,s,ance from the transmitter; applying an electromagnetic (EM) wave field to 
>hc strata using the transmitter; delecting the EM wave field response using me 
rece.ver; extracting phase infoimafion from me wave response; repea.ing .he 
procedure with the transmitter and/or receiver in different localions for a plurality of 
.ransm,ss,o„ S ; and using «h e phase information from «he wave response for the 
plurahly of transmissions, in order to determine the presence and/or nature of any 
reservo.r present; in which melhod, the wavelenglh of me transmission is given by the 
formula 



0.1s < X< 10s; 



wherein X is the wavelength of the transmission Ihrough Ihe overburden and s is the 
distance from Ihe seabed ,o the reservoir; the predetermined offset between the 
transmitter and a receiver is given by the formula 



< L < 10X; 



where X is the wavelength of the transmission through the overburden and L is the 
distance between the transmitter and the receiver; and the transmission frequency is 
from 0.01 Hz to 1 kHz. 



Preferably, the wavelength of the transmission is in the r 



ange 



0.1s<A.<5s. 



More preferably A is from about 0.5s to 2s. The transmission fluency is preferably 



from 0. 1 to 20 Hz, for example I Hz. 

Thus the offset can be varied by moving ,he receiver; or tadee( , ^ ^ o[ 
even oft. Altemativeiy, the pre4e,ennmed offset can be kept constant by m ovi„ g 
^ both the transmitter and receiver. 

consent representing a refracted wave, and then extracting and ustng the refracted 
wave response. 
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Thus, the honzonta, boundaries of the reservoir may be found by analysing to sIope 
and/or slope change of the cnrvefs, of phase and/or magnitude as a function of SO urce- 
recerver offset distance or position, or by ana,ysing the variation in phase and/or 
ma^ttude for a frxe, source-receiver offse, a, sevem, Nation, The most useful 
source.rece.ver offse, is typically !ar g er tban the -cnhcal offset". Fn this par, of the 
curve, the change in slope, may indicate the reservoir boundary. 

Both the source and the receiver are preferably inside the reservoir area to achieve the 
ma„es slope (or ^adient). This is true for both the phase and the magnitude curves 

20 Soon afrer e,ther the source or the receiver ,eaves the reservoir area, the shape 
-crease rap,d.y. From the position where this change occurs, the reservoir bounly 
may he mapped. The true reservo.r boundary wit, probably he closer the centre of Z 
— compared to the .ocation where the s.ope change occurred, typically ,0 to 
20/o of the reservoir depth. The detaited position may be caicu.ated using the 

25 measured data and forward modelling. 

This technique of the invention can be used in conjunction with conventional seismic 
teohmques to ■dentify hydrocarbon reservoirs. 

30 If the offse, between the transmitter and receiver is significantly greater than three 



mes the depth of ,hc reservoir f rom the se abeo (i.e. the thickness of the overburden) 
> t «D he appreciated tha, the attenuation of the refracted wave wiH often he iess than' 
of direct wave and the reflected wave. The reason for this is the fact tha, the p ath 
of the refracted wave win he effectively distance from the transmitter down , 0 the 
reservotr t.e. the thickness of the overburden, p,us the offset aiong the reservoir pl 
the .stance front the reservoir up ,„ the receivers u once again fc to . 



overburden. 



us 
the 



f no H/C reservoir is present in the area of the transmitter and receiver, the phase of 
he detected wave response win consist of a direct wave and wil, therefore change 

hneafly with a changing offset. S.mi.arly, the phase of the detected wave response 

wtll reman, constant a, a constant offset in different locations. 

However, if an H/C reservoir is present, mere wiU be a retracted wave component in 

the wave response and this may prrfommate. Due to the higher phase vetocity 

(wavespeed) in H/C fined strata, this wih have an effect on the phase of the received 

wave response. .„ the case of an increasing offset, the phase wifl no, change linearly; 

a Plot o phase against offset wi„ be a curve with a constanfly changing s!ope. Thus 
change ft„ ra . straight Une t0 , ^ ^ ^ ^ ^ ^ 

H/C reservoir. 



a 
an 



In the case of a constant offset, the presence of an H/C rese™ir wil, give rise to a 
constant hu, different phase va.ue at different locations compared ,o ,he situation 
where no H/C resen-oir is present. Thus, a change m phase value wil. indica.e the 
boundary of an H/C reservoir. 

Preferably, (he procedure is repealed at different offsets. 

m one embodiment, the method mc.udes p.otting a graph of the phase of the refracted 
wave response from a particular stratum against offae, and amdysmg the slope of the 
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graph in order ,c determine the nature of the stratum. AltemativeIy , , he ^ 

* P 0ttme ' graP " ° f " K ^ ° { «* -fl-ed wave response front a parocu.ar 
stratum and .dentifying a ehange in the slope of the graph. 

In another embodiment, the method comprises repeating the procedure a, different 
ocatmns, using the same predefined offset, and analysing the phase of the 

z?ir response from a partia " ar straium in t,,der *° ^ » *■» *> 

The poiarization of the source transmisston win determine how much energy is 

irr T T th : oi '' beann8 ,ayer ,n the direc,ion ° f ,,k a - w 

s therefore the seiected transmitter. ,„ genera,, i, „ pr eferab,e to adopt a dipole with a 
arge effeefve .ength. The transmitter dipoie may therefore be , 00 to 1000 meters in 
engdr and may be towed in two orthogona. directions. The receiver dipo,e optimum 
length is determined by the thickness of the overburden. 

The techniuue can be apphcab.e in explonng ,and-based subterranean reservoirs bu, is 
espectady applicable to snbmarine, in particuiar sub-sea, subterranean reservoirs 
Preferable field , appfied using one or more transmit ,„ cate d on me earth's 
urfece, and the detection is carried on, by one or more receivers iocated on the ea^s 

ctrr , a p r md appiication ' ,he - on or 

close to the seabed or the bed of some other area of water. 

The transmitted f,e,d may be puiscd, however, a coherent continuous wave optional* 
wtth stepped fluencies is preferred. „ may be transmitted for a si^ifican, period of 
tune, durtng which .he transmutes should preferably be stationary (although i« cou.d 
be movmg sI ow ly) , and the transmission stabie. Thus, the f,e,d may be transmitted for 
a penod of tnne from 3 seconds ,„ 60 minutes, preferably from , 0 seconds ,o 5 
mmu.es, for example about , minute. The receivers may also be arranged to detect a 
dtrec, wave as well as the refracted wave from the reservoir, and the ana!ysis may 
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15 



20 



include extracting phase and amptitude data of the refracted wave from corresponding 
data from the direct wave. 

I. wil! be appreciated that the present invention may be used to de.en.ine the position 
he extent, ,he nature and the volume of a particniar stmtum, and may a,so be used to 
detect changes in these parameters over a period of time. 



measures 



25 



30 



The present invention also extends to a method of surveying subterranean 
winch comprises; performing a seismic survey to determine the geological sZZ of 
a regron; and where that survey reveals the presence of a subterranean reservoir 
subsequently performing a method as described above. 

The invention may be carried into practice in various ways and win now be iHustrated 
tn the followmg simulated Examples. In the accompanying drawings. 

Figure l is a schematic representation of a situation in which measurements are taken 
outside the area of a reservoir; 

Figure 2 is a schematic representation of a situation in which measurements are taken 
mside the area of a HC reservoir; 

Figure 3 is a schematic representation of a situation in which measurements are taken ' 
crossing a H/C reservoir boundary; 

Figure 4 is a view similar to Figure 3, in which the receiver is in a borehole; and 
Figure 5 is a graphiea! representation of the results from the situations in Figures 1 to 



In Figures 1 ,o 4, i, is assumed that the sea III—, resistivity of 0.3 ohm-m and a 
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depth of 1 000m from the surface 1 2 to the seah^H n n 

Ine sea oed 13. The overburden 14 has a 

resrstivity of 0.7 ohm-m. The transmitted wave had a frequency of 1 Hz. 

Figure I shows the station where a receiver 1 5 and a transmitter 16 are located on 
the seabed 13 in a region where there is no subterranean reservoir The position of the 
receiver IS. fixed. The transmitter ,6 „ moved to various different positions 
resulting in different offsets from the receiver 15. At each position of the transmitter 
, an EM field is transm.tted and the wave fie.d response is detected by the receiver 
5. The wave field response mclndes a direct wave component and reflected wave 
component, Phase information is extracted for each response and the results are 
Plotted as a graph in Figure 5. The p,o, of phase against offset (source-receiver 
stance) is a straight line 51, showing tha, the phase varies linearly with offset. 

Figure 2 shows the situation where the transmitter 1 5 and receiver ,6 are depmyed in a 
regton where mere is a subterranean H/C reservoir . 7. The reservoir , 7 is , 00m thick 
and has a resistivity of 50 ohm-m, and is located a, a depth of 800 m beneath the 
seabed 13. The procedure described in relation to F.gure 1 is repeated and again the 
results are shown in Figure 5. 

In this case, the wave field response additional mdudes a refracted wave component 
from the reservoir ,7. This affects the phase of the response and resu,« S in the plo, of 
Phase against offset talcing the form of a smooth curve 52, with a constant change of 
slope. ° 



Ftgure 3 shows the situation where tire transmitter 15 and receiver 16 are deptoyed in a 
regton where there is an H/C reservoir , 8 with a boundary 19. The reservoir 1 8 is 
suntiar to that in F.gure 2 bu, its boundary 19 is a, a position 3km away from the 
portion of the transmitter 1 5. Again, the procedure described in re.ation ,o Figure 1 is 
repeated and the results are shown in Figure 5. 
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In On, case, ,he wave field response phase curve 53 initially follows toe curve 52 of 
•he results from Figure 2 bu, soon after an offset of 3km, the results follow a straight 
ime whtch has the same slope as the line 5 1 of the Figure 1 situation. The change in 
slope of the curve 52 clearly shows the presence of a boundary between an H/C 
reserve, and no H/C reservoir. The posit.on of the boundary ,9 can be deduced from 
the position of the change in slope. 

The situation shown in Figure 4 is sitni.ar to tha, of Figure 3 except that in this ease 
fte receiver 25 is located in a we,, bore 27 and the transmitter 26 is moved to differ™, 
ocattons to vary the offset. Once again, the procedure described in relation to Figure 
1 was repeated and the results ate shown in Figure 5. 

Due to the location of the receiver 25 in a well bore 27 and within the reservotr 18 the 
Phase curve 54, as the offset increases, follows a straight line initiaUy. The slope of 
the mtttal portion is shallower than the slope of the curve 5 1 where no H/C reservoir 
was present, due to the effect of the presence of a direct wave propagated through the 
H/C reservoir. However, soonafter an offset of 3km, the slope of the curve 54 begins 
«o change and continues ,o change until i, attains the same slope as the curve 5 1 of 
Ftgure , . This change again shows clearly the presence of an H/C reservoir boundary 
and tts position can be deduced from the position of the change in the slope of the 
curve 54. 
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Claims 



1. A method of investigating subterranean strata which comprises: deploying an 

electric dipole transmitter antenna; deploying an electric dipole receiver antenna at a 

predetermined offset distance from the transmitter; applying an electromagnetic (EM) 

wave field to the strata using the transmitter; detecting the EM wave field response 

using the receiver; extracting phase information from the wave response; repeating the 

procedure with the transmitter and/or receiver in different locations for a plurality of 

transmissions; and using the phase information from the wave response for the 

plurality of transmissions, in order to determine the presence and/or nature of any 

reservoir present; in which method, the wavelength of the transmission is given by the 
formula 



0.1s < X < 10s; 



wherein A is the wavelength of the transmission through the overburden and s is the 
distance from the seabed to the reservoir; the predetermined offset between the 
transmitter and a receiver is given by the formula 



0.5k <L< 10X; 



where X is the wavelength of the transmission through the overburden and L is the 
distance between the transmitter and the receiver; and the transmission frequency is 



from 0.01 Hz to 1 kHz. 



2. 



A method as claimed in Claim 1 , in which the procedure is repeated at different 
offsets. 



3. A method as claimed in Claim 2, which includes plotting a graph of the phase 
of the reflected wave response from a particular stratum and identifying a change in 



the slope of the graph 
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4. A method as claimed in Claim 1 , which comprises repeating the procedure a, 
Afferent lo ca,ions, using „ e same predetennined ^ ^ fe rf 

any refracted wave response from a parti cu,ar stratum in order to identify a change in 
the phase value 

5. A method as claimed in any preceding Claim, which additionally includes 
extmctmg and using amplitude information from the wave response. 

0- A method as claimed in any preceding Claim, in which the transmitter and/or 
recerver ,s located on or close to the seabed or the bed of some other area of water. 

7 - A method as claimed in anv orecedina rioi™ • u t , „ 

. any P recedl ng Claim, m which the frequency of the 

EM field „ contmuously varied over the transmission period. 

* A method as claimed in any preceding Claim, in which the field is transmitted 
for a period of time from 3 seconds to 60 minutes. 

9- A method as claimed in Claim 8, in which the transmission time is from 10 
seconds to 5 minutes. 

fl A T d " C ' aimed PreCedin8 C,am ' » "** «» — sio„ 

frequency is from 0. 1 to 20 Hz. 

11. A method as claimed in any preceding Claim, which includes suppressing any 
drrect wave and/or any other ta wn wave contribution that may disturb th 
measurements, thereby reducing the retired dynamic range of the receiver and 
increasing the resolution of any refracted wave. 
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12. A method of surveying subterranean strata which comprises: performing a 
seismic survey to determine the geological structure of a region and where that survey 
reveals the presence of a subterranean reservoir, subsequently performing a method as 
claimed in any preceding Claim. 

1 3. A method of producing a survey report which comprises carrying out a method 
as claimed in any preceding Claim and generating a report based on the results of the 
investigations. 
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